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Abstract. We present a nucleosynthesis sensitivity study for the y-process in a Supernova type II model within 
the NuGrid research platform. The simulations aimed at identifying the relevant local production and destruc¬ 
tion rates for the p-nuclei of molybdenum and at determining the sensitivity of the final abundances to these 
rates. We show that local destruction rates strongly determine the abundance of ®^Mo and ^''Mo, and quantify 
the impact. 


1 Introduction 

The p-nuclei between ^'^Se and can be produced un¬ 
der explosive conditions in a sequence of photodissocia¬ 
tion reactions and subsequent y6-decays starting at s- and 
r-process seeds |[ll]. Most of the p-nuclei are about two 
orders of magnitude less abundant than other stable iso¬ 
topes belonging to the same element. Relevant excep¬ 
tions are the isotopes and Their pro¬ 

duction in stars is a puzzle for nuclear astrophysics since 
present models underproduce these p-nuclei by orders of 
magnitude. According to recent stellar model calculations, 
®"^Mo is mainly synthesized via the (y.n) photodisintegra¬ 
tion chain starting from the more neutron-rich and stable 
molybdenum isotopes Some of the light p-nuclei, like 
the neutron magic isotope ®^Mo, may also be synthesized 
by proton capture reactions |0l- 

We investigate the nucleosynthesis of ®^Mo and ®^Mo 
in a Supernova type II model. We present results on the 
sensitivity of rates of direct production or destruction re¬ 
actions of the two Mo isotopes. 

2 PPN simulations using a Supernova 
type II model 

For the post-processing nucleosynthesis (PPN) simula¬ 
tions we use a stellar progenitor with an initial mass of 
25 Mq and metallicity Z = 0.02 ||2l. The p-process nucle¬ 
osynthesis is calculated by using classic Supernova type II 
trajectories 0. In the model, a shock front passes through 
the Ne/O burning zone, which is subdivided into 14 mass 
layers. Figure [T] shows the temperature profiles, which 
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describe the astrophysical environment of each layer as 
a function of time. The innermost mass shell is the 
hottest and densest environment with a peak temperature 
of 3.45 GK. The maximum value drops to 1.79 GK for the 
outermost mass layer. The shock front reaches the mass 
layers successively. The temperature rapidly increases to 
the maximum value and afterwards drops slowly. 
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Figure 1. Temperature profiles of the classic Supernova type II 
trajectories @]. Left to right: profiles from the innermost mass 
layer to the outermost. 


The complete nucleosynthesis is calculated using the 
post-processing code PPN within the NuGrid research 
framework ITllst]. NuGrid offers a software framework 
for nucleosynthesis simulations in relevant astrophysical 
environments. A large network of over 5,000 isotopes and 
more than 60,000 reactions is used. Most of the nuclear 
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Figure 2. Time-integrated nucleosynthesis fluxes producing and 
destroying the isotope ^^Mo. The sum of the production fluxes is 
normalized to 100%, and the destruction fluxes are scaled with 
the same factor. Fluxes smaller than 1% are not shown. 


reaction rates adopted are taken from the JINA Reaclib 
Database Vl.l lot]. 

3 Nucleosynthesis fluxes 

The relative time-integrated nucleosynthesis fluxes of all 
mass layers for the isotopes ®^Mo and ^'^Mo are shown in 
Figures |2] and |3] For each isotope, the sum of the pro¬ 
duction fluxes is normalized to 100%, and the destruction 
fluxes are scaled with the same factor. Fluxes smaller than 
1% are not shown. If the sum of destruction fluxes is less 
than 100%, there is a net yield in the simulation. Other¬ 
wise the isotope gets depleted. The sum of the production 
fluxes is larger than the sum of the destruction fluxes for 
both isotopes. The isotope ®^Mo has the largest net yield 
with 100% production and 93% destruction. 

The isotope ®^Mo is mainly produced by photodisin¬ 
tegration of ^^Mo. The reaction ®^Mo(y,p) is the main 
destruction path followed by ^^Mo(n,y) reactions. The 
reaction ®^Mo('y,n) strongly contributes to the production 
fluxes of ^"^Mo. The isotope ®^Mo is equally destroyed by 
( 7 ,n) and (n,y) reactions. 

The reactions accounting for a significant fraction of 
the nucleosynthesis fluxes are investigated in the rate vari¬ 
ation study in section |4] 

4 Abundances of and 
depending on local production and 
destruction rates 

The PPN simulations using a Supernova type II model aim 
at identifying the relevant reactions for the production and 
destruction of the Mo p-nuclei, and at quantifying the ef¬ 
fect of rate changes on the final abundances. Previous nu¬ 
clear sensitivity studies varied the reaction rates for all re¬ 
actions of a certain type in the network, as well as for all 
reactions involving the p-nuclei simultaneously fill. In 
this work, we explore the impact of single rates directly 
producing or destroying ®^Mo and ^"^Mo. 

The effect of single rate variations on the final abun¬ 
dances of the p-nuclei ®^Mo and ®^Mo was determined 
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Figure 3. Same as Figure|2l for 


by multiple simulations of the same scenario using differ¬ 
ent reaction rates. The rates were multiplied/divided by 
factors 2 and 5 for reactions involving a neutron ((y,n) 
and (n,y)) or a proton ((y,p) and (p,y)), and by factors 
5 and 10 for reactions involving an or-particle ((■Y,a) and 
The factors were chosen according to the com¬ 
monly stated uncertainties of the reaction rates (fi] and 
references therein). In the PPN simulation, the rate taken 
from the library is calculated for the current temperature in 
the astrophysical environment and then multiplied by the 
factor set by the user. 

This work shows the dependencies for the rates where 
the maximum abundance ratio is larger than 5% within the 
applied rate variations. 

4.1 The case of ®^Mo 

Figure |4] shows the final abundance of ®^Mo as a func¬ 
tion of the factor applied to one production or destruc¬ 
tion rate. The abundance of ®^Mo strongly depends on 
the ®^Mo(y,p) reaction rate, which is the main destruction 
path. The abundance of ®^Mo drops with increasing rate 
and vice versa (Figure|4] top panel). Neutron dissociation 
reactions are the main production path of ®^Mo (Figure |4] 
bottom panel). The ^‘^Mo(y,n) reaction shows a strong im¬ 
pact on the abundance of ®^Mo. Variations of the reaction 
rate of ®^Mo(y,n) hardly affect the ®^Mo abundance since 
the unstable isotope ^^Mo has to be produced first. 

The ®^Mo abundance drops if the reaction rate of 
®'*Mo(y,Qf) increases. The ratio of the reaction rates of 
®"^Mo(y,a) and ®"^Mo(y,n) is larger at high temperatures, 
hence, the nucleosynthesis path is altered producing less 
®^Mo via (y,n) reactions. 

4.2 The case of ®^l\/lo 

Photodisintegration reactions starting from heavier Mo 
isotopes produce ®^Mo. The destruction reaction 
®"^Mo(y,n) mainly determines the final abundance of ®^Mo. 
In the simulation, a small fraction of the nucleosynthe¬ 
sis fluxes proceeds via ®'*Mo(y,a). An increase of the 
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Figure 4. Impact on the abundance of ®Mo depending on the 
y-induced production and destruction rates. The results for the 
reaction rate of ^^Mo( 7 ,p) (top) are shown in a single plot due to 
its large impact on the abundance compared to the other reaction 
rates. 



Figure 5. Impact on the abundance of ^'^Mo depending on the 
y-induced production and destruction rates. 


®'*Mo(y,a) rate reduces the abundance of signifi¬ 

cantly. If the ®^Mo(y,n) rate is increased, more ®^Mo is 
destroyed and not available for the production of ^"^Mo via 
®3Mo(n,y). 


5 Summary and outlook 

The impact of variations of local production and destruc¬ 
tion rates on the p-nuclei ®^Mo and ^"^Mo was investi¬ 
gated. The abundance of ®^Mo is mainly determined by 
the ®^Mo(y,p) reaction rate, the ®"^Mo abundance by the 
^"^Mo(y,n) and ®‘*Mo(y,a) reaction rates. Most rates taken 
from JINA Reaclib Database V1.1 are calculated within 
the Hauser-Feshbach model New evaluated data can 
be compared to these rates and new final abundances of 
the Mo p-nuclei can be obtained from the results presented 
here. 

In the y-process, reactions on heavy nuclei (Z > 42) 
may determine the reaction flow to the light p-nuclei and, 
hence, the abundance of ®^Mo and ®"^Mo. Future studies 
foresee a global sensitivity study of reaction rates in dif¬ 
ferent p-nuclei production scenarios including Supernovae 
type la. 
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